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The effects of adiabatic phonons on a spin-fermion model for high T c cuprates are studied using 
numerical simulations. In the absence of electron-phonon interactions (EPI), stripes in the ground 
state are observed [1] for certain dopings while homogeneous states are stabilized in other regions of 
parameter space. Different modes of adiabatic phonons are added to the Hamiltonian: breathing, 
shear and half-breathing modes. Diagonal and off-diagonal electron-phonon couplings are consid- 
ered. It is observed that strong diagonal EPI generate stripes in previously homogeneous states, 
while in striped ground states an increase in the diagonal couplings tends to stabilize the stripes, 
inducing a gap in the density of states (DOS) and rendering the ground state insulating. The off- 
diagonal terms, on the other hand, destabilize the stripes creating inhomogeneous ground states 
with a pseudogap at the chemical potential in the DOS. The breathing mode stabilizes static di- 
agonal stripes; while the half-breathing (shear) modes stabilize dynamical (localized) vertical and 
horizontal stripes. The EPI induces decoherence of the quasi-particle peaks in the spectral functions. 

PACS numbers: 74.20.De, 74.25.Kc, 74.81.-g 
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I. INTRODUCTION 

The pairing mechanism responsible for high T c super- 
conductivity is still unknown. The electron-phonon inter- 
actions that satisfactorily explain pairing for traditional 
superconductors within the BCS theory [2] would require 
phonon frequencies incompatible with the material sta- 
bility in order to produce the observed high critical tem- 
peratures in the cuprates. [2] For this reason many re- 
searchers believe that magnetic interactions, which are 
observed in all the cuprates, may play an important role 
in the pairing mechanism. [3] As a result of this hypoth- 
esis, most of the Hamiltonians proposed to study the 
physics of the cuprates, such as the Hubbard and t-J 
models, only incorporate electronic and magnetic degrees 
of freedom. [4] However, experiments indicate that there 
are active phonon modes in the cuprates. [5-8] In addi- 
tion, experimentally a very rich phase diagram, particu- 
larly in the underdoped regime, is started to be unveiled. 
Some materials appear to have ground states with elec- 
tronic stripes as observed in neutron scattering results 
[8], while scanning tunneling microscopy indicates nano- 
size patches of superconducting and non-superconducting 
phases. [9,10] The emerging phase complexity is remi- 
niscent of the experimental data for manganites where 
competing electronic, magnetic and phononic degrees of 
freedom are responsible for the rich phase structure. [11] 

For these reasons it is important to include electron- 
phonon interactions in models for the cuprates. This 
would allow to understand whether EPI stabilize or 
destabilize charge stripes, what kind of inhomogeneous 
textures, if any, develop and, eventually, whether the 



interplay of magnetic and phonon interactions with the 
electrons is responsible for the pairing mechanism. [12] 

The first step towards the goal of introducing EPI in 
models for the cuprates is to propose a simple but physi- 
cally realistic Hamiltonian that can be studied with unbi- 
ased techniques. The proposals already in the literature 
include momentum dependent electron-phonon couplings 
that lead to long-range interactions in coordinate space 
[13,14] and/or quantum phonons which are very difficult 
to treat numerically. [15,16] Most studies have been per- 
formed using mean-field, slave-boson, or LDA approxi- 
mations. [17,18] Numerical simulations have been done 
using the t — J model in very small lattices, with a lim- 
ited number of phonon modes and diagonal couplings, 
[15,16] or on the one-dimensional Hubbard model. [19] 

In this paper we will study numerically a spin-fermion 
(SF) Hamiltonian for the cuprates [1] with electron- 
phonon interactions. This model reproduces many prop- 
erties of the cuprates and presents stripes in the ground 
state, due solely to spin-charge interactions, in some re- 
gions of parameter space. [1] Thus, it provides a frame- 
work particularly suitable to study the effects of electron- 
phonon interactions on the preformed stripes. However, 
charge homogeneous ground states are also found in other 
regions of parameter space which allows to investigate 
charge inhomogeneity induced by EPI. Several phononic 
modes will be studied and diagonal and non-diagonal 
couplings, i.e., the dependence of the hopping and other 
Hamiltonian parameters on the lattice distortions, will be 
considered. The work will be performed in the adiabatic 
limit, i.e., at zero phononic frequency. 

The paper is organized as follows: in Section II the 
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Hamiltonian is introduced; the effect of electron-phonon 
interactions on striped states are presented in Section 
III, while Section IV is devoted to the effects of EPI on 
homogeneous states. Section V contains the Conclusions. 

II. THE MODEL 

The SF-model is constructed as an interacting system 
of electrons and spins, mimicking phenomenologically the 
coexistence of charge and spin degrees of freedom in the 
cuprates [20]. Its Hamiltonian is given by 

H = -* E ( c L c j« + h - c -) + J E s; • Si + J ' E S; • S J' 

<ij>" 1 <ij> 

(1) 

where C; Q creates an electron at site i = (i x , i y ) with spin 

projection a, Si=^2 a p c\ a &a0 c ip 1S the s P m of the mo- 
bile electron, the Pauli matrices are denoted by a, Si is 
the localized spin at site i, (ij) denotes nearest-neighbor 
(NN) lattice sites, t is the NN-hopping amplitude for the 
electrons, J > is an antiferromagnetic (AF) coupling 
between the spins of the mobile and localized degrees of 
freedom, and J' > is a direct AF coupling between the 
localized spins. The density (n) = l — x of itinerant elec- 
trons is controlled by a chemical potential fi. Hereafter 
t = 1 will be used as the unit of energy. J' and J are fixed 
to 0.05 and 2.0 respectively, values shown to be realistic 
in previous investigations [1]. The temperature will be 
fixed to a low value: T=0.01, which was shown before to 
lead to the correct high-T c phenomenology. [1,21] 

The diagonal electron-phonon part of the Hamiltonian 
being proposed here is given by 

^e ( i ) P h = - A E^ ) ^' ( 2 ) 

i 

where ri\ — ^ CT c\ a c\ a is the electronic density on site 

i and Q\ ' is the phonon mode defined in terms of the 
lattice distortions u- ha which measures the displacement 
along the directions a = x or y of oxygen ions located at 
the center of the lattice's links in the equilibrium posi- 
tion, i.e., u- ha — 0. The index (j) identifies the phonon 
mode. In this work, the following phonon modes will be 
considered: 

(a) The breathing mode given by 

<9i (1) = E( Ui < a - u i-&,a); ( 3 ) 

a 

(b) The shear mode, in which the oxygens along e.g., x 
move in counterphase with the oxygens along y, given by 

Ql (2) =£(-l) <T («l,a-Ul-a,a), (4) 

a 

with a = 1(— 1) for a. = x(y); 



(c) The half-breathing mode along x given by 

Qf ] = {ui,x - Mi-x,x); (5) 
and (d) The half-breathing mode along y given by 

<3i (4) = («!,„ -«!-*,»)■ (6) 

Note that although the proposed interactions seem 
local in coordinate space, they correspond to cooper- 
ative lattice distortions which, in turn, will produce 
strongly momentum dependent effective electron-phonon 
couplings, in agreement with the experimental evidence 
observed in the cuprates. [6] 

A term to incorporate the stiffness of the Cu-0 bonds 
is added. The term bounds the amplitude of the lattice 
distortions induced by i/ c „ p h. Its explicit form is: 

H ph = K^2(ui. a ) 2 , (7) 

i,a 

where k is the stiffness parameter that will be set to 
1 here. In addition, we will consider the off-diagonal 
interactions induced by the lattice distortions. To ob- 
tain these terms we follow the approach of Ishihara et al. 
[19,18] As a result the hopping t in Eq.(l) now becomes 
site and direction dependent and it is given by 

fij=* + 7[«(i)+u(j)]. (8) 
where 7 is a parameter and 

The Hciscnbcrg coupling J' in Eq.(l) also is affected by 
the lattice distortions and it has to be replaced by 

Jii = J' + 9.n[u(i) + u(3)}, (10) 

where gj is another parameter. 

As stated above, the spin-fermion model with electron- 
phonon interactions will be studied with a Monte Carlo 
(MC) algorithm. To simplify the numerical calculations, 
avoiding the sign problem, the localized spins are as- 
sumed to be classical (with |Si|=l). This approximation 
is not drastic since most of the high Tc phenomenology 
is reproduced in this limit, and it was already discussed 
in detail in Ref. [1]. Details of the MC method can be 
found in Ref. [22]. Square lattices with 8x8 and 12 x 12 
sites will be studied here. 



III. INFLUENCE OF PHONONS ON STRIPED 
STATES 

Neutron scattering experiments have shown that dop- 
ing causes magnetic incommensuration in the cuprates. 
[23,24] The origin of this phenomenon is still being de- 
bated. One possible scenario is the formation of charge 
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stripes in the ground state upon doping. [25,26] Experi- 
ments on nickelates such as La2- x Sr x Ni04 (LSNO) have 
shown the presence of diagonal static stripes [8] and there 
is evidence of stripes also in LNSCO, i.e., Nd doped 
La 2 „ x Sr x Cu0 4 (LSCO) [27] and in La 2 _ x Ba x Cu0 4 
(LBCO). [28] It is conjectured that the magnetic incom- 
mensurability observed in other high T c cuprates is due to 
the presence of dynamical stripes and that the stripe dy- 
namics may be related to the electron-phonon couplings 
in the different materials. Despite these theoretical sce- 
narios, it has been very difficult to find striped ground 
states in models for the cuprates when studied with unbi- 
ased techniques. Stripes have been observed in the t — J 
model [29,30] but they are difficult to stabilize and it is 
not clear whether the striped state is the actual ground 
state or a very low lying excited one. Stripes, on the other 
hand, have been obtained without biases in the SF model 
that will be studied in this work. [1] This characteristic 
will allow us to explore the effect of electron-phonon in- 
teractions in striped ground states. The charge texture 
will be monitored by measuring the density structure fac- 
tor N(q). [1] The magnetic structure factor S'(q) will 
provide information on the magnetic properties. Lattice 
distortions will be monitored by measuring correlations 
between the displacements, although the development of 
long range correlations in the lattice degrees of freedom 
is not anticipated. 




FIG. 1. (a) MC snapshot of an 8 x 8 lattice at (n) =0.875 
for A = and 7 = 0. The size of the circles is proportional to 
the electronic density; the shaded circles have charge density 
larger than the average, i.e., m > (n) = 0.875. The arrows 
represent the projection of the localized spins in the plane 
x — y; the lines indicate lattice distortions (see text); (b) same 
as (a) but for A = 1 and mode Q^; (c) same as (b) but for 
A = 2. 

The dynamic properties of the system will be moni- 
tored by measuring the density of states N(uj), the one- 
particle spectral functions A(q, uS) and the optical con- 



ductivity <j(ui). [1,21] 

A. Diagonal electron-phonon term: 

In our investigations it has been observed that in gen- 
eral the diagonal electron-phonon interaction plays an 
stabilizing role on stripe structures. This behavior was 
obtained for the four phonon modes studied here. For 
< A < 2 the holes become more localized in the stripes 
as A increases. This can be seen in Fig.l, where snap- 
shots for (n) = 0.875 are displayed for A = (Fig. la), 
A = 1 (Fig. lb), and A = 2 (Fig.lc). The lines in the 
snapshots indicate the lattice distortions. If all the dis- 
placements Ui jQ were then, the lines would cross at the 
middle point of the links that join the lattice sites (as in 
Fig. la). Out of center crossings indicate ionic displace- 
ments. It is clear from the figure that, as A increases, 
lattice distortions in the direction perpendicular to the 
stripe develop along the stripe with the mode Q^ 2 \ fur- 
ther localizing it. In Fig.lc large displacements along the 
horizontal direction can be seen in the links next to the 
stripe. The stripes become thinner and the density of 
holes per site inside the stripes increases. 




FIG. 2. (a) MC snapshot of an 8 x 8 lattice at (n) =0.875, 
for mode A = 2 and 7 = 0; (b) same as (a) but for mode 

Q^ 2 '; (c) same as (a) but for mode Q' 3 '. 



We have observed some differences between the effects 
of the various phonon modes studied here as the strength 
of the diagonal electron-phonon coupling A increases. In 
Fig. 2 snapshots for A = 2 at (n) — 0.875 are presented for 
Q^ 2 \ and Q( 3 \ [31] A clear tendency to form diag- 
onal stripes is seen for the breathing mode (Fig. 2a). 
Since in this case holes are localized by being surrounded 
by four elongated bonds they cannot be accommodated 
in vertical or horizontal formations. Note that the single 
stripe that is stable at A = for the electronic density 
shown in Fig. 2 gets destabilized due to the QW strong 
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electron-phonon coupling. This result agrees with the 
fact that diagonal stripes are observed in LSNO and ex- 
periments indicate that the breathing mode is the mode 
most strongly coupled to the electrons. [32] According to 
our results, a robust diagonal coupling of the electrons to 
the breathing mode should be expected in the nickelates. 

In Fig. 2b it can be observed that the shear mode 
tends to stabilize vertical (or horizontal) stripes because 
a large horizontal (or vertical) distortion occurs, localiz- 
ing the holes along the stripe. Interestingly, the half- 
breathing mode also produces vertical (or horizontal) 
stripes but the holes are less localized since the lattice can 
distort only along the horizontal (or vertical) direction. 
As a result, more dynamical stripes are observed for the 
half-breathing modes even for strong diagonal electron- 
phonon couplings (Fig. 2c). Notice that the experimental 
evidence indicates that in LNSCO, where vertical and 
horizontal stripes are observed, the mode more strongly 
coupled to the electrons is the half-breathing mode. [6] 

Another indication of increasing localization with in- 
creasing A is observed in the peak in N(q), at q = 
("7r/4, 0) for the stripes shown in Fig.la-c, which becomes 
better developed as A increases (Fig. 3a). In addition, the 
pseudogap in the density of states at the chemical poten- 
tial becomes a full gap indicating an increase in insulating 
behavior (Fig. 3b). 



(b) 




N(co) 



(n,0) 

q 



(0,0) 




-0.9 



FIG. 3. (a) The charge structure factor iV(q) for various 
values of A for the same parameters as in Fig.l; (b) The 
density of states N(u>) for several values of the diagonal elec- 
tron-phonon coupling A, for the same parameters as in (a) The 
phonon mode is . 

Similar results are obtained at (n) — 0.75 for which 
two stripes are stabilized, even at A = 0, in the 8x8 
systems studied here. 

The general trend, for different electronic densities and 
phonon modes, is that the ground state becomes more 
insulating as the diagonal electron-phonon coupling in- 
creases. The shear mode will be used as an example 
but similar qualitative behavior is observed for the other 
modes. In Fig. 4a it can be seen that the spectral weight 
at oj — in the density of states decreases with increas- 
ing A for different values of the electronic density (n). 
The Drude weight, shown in Fig. 4b, also decreases and 
insulating behavior is obtained for all densities at A = 2. 

Although the effect of charge localization with increas- 
ing A is more pronounced at densities for which stripes 



are observed, we see in the snapshots shown in Fig. 5 for 
(n) = 0.8 and mode that the charge becomes more 
localized as A increases and the lattice distortions localiz- 
ing the holes develop large values for A w 2. Note that AF 
domains separated by walls of holes are observed. These 
inhomogeneous structures appear to replace the stripes 
when the density of holes is not commensurate with the 
lattice size. An important characteristic of these inho- 
mogeneous states is that, although no features are ob- 
served in the charge structure factor AT(q), incommensu- 
rate magnetic correlations are still present and the peaks 
in S(q) occur at (tt, tt — 5) and (tt — S, tt), i.e., in qualita- 
tive agreement with the incommensurate peaks observed 
in the cuprates. 
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FIG. 4. (a) Spectral weight in the density of states N(u) at 
ui — as a function of the diagonal electron-phonon coupling 
A for several values of the electronic density (n) and mode 
Q' 2 '; (b) Drude weight as a function of the diagonal elec- 
tron-phonon coupling A for several values of the electronic 
density (n) and mode Q^ 2 ' . 

In the configurations shown in Fig. 5, S(q) has a maxi- 
mum at q = (tt, 37r/4) for snapshot (a), but there is also 
a less intense peak at q = (37r/4, tt). 

Incommensurate peaks in S(q) at q = (tt, 37r/4) and 
q = (37r/4,7r) with almost equal weight are observed for 
(n) = 0.8 and A = 1 (see Fig.5d) although no peak in 
N(q) is observed. In Fig. 5b the snapshot of the final 
configuration of the corresponding Monte Carlo run is 
shown. Other snapshots of configurations appearing dur- 
ing the measuring part of our simulation are displayed 
in Fig. 6. It can be seen that the ground state is not 
frozen and that there is a dynamical charge redistribu- 
tion. Thus, the magnetic incommensurability observed 
in some cuprates could be due to charge and spin config- 
urations similar to those presented in Fig. 6, which could 
be interpreted as "dynamic stripes" . In fact, these results 
may indicate that the patches in Fig. 5 and 6 are not ran- 
dom islands since, if that were the case, we would expect 
that the maxima in S(q) would form a ring in momen- 
tum space. [33] It is tempting to associate the observed 
patches with "dynamical" stripes. Notice that for the 
states with "static" stripes only one peak at q = (tt, tt—S) 
is observed in 5(q) if the peak in N(q) is at q = (0, 2(5). 
On the other hand, the states with "dynamical" stripes 
naturally reproduce the four peaks observed in neutron 
scattering experiments for the cuprates [23,24] and the 
"patch-like" shape of the clusters would explain, at the 



4 



same time, the apparently random inhomogeneous struc- 
tures observed in STM experiments. [9,10] 

As the diagonal electron-phonon coupling A increases 
beyond A = 2, important quantitative changes are ob- 
served, in particular, the tendency to bipolaron forma- 
tion induced by strong lattice distortions. This behavior 
will be discussed in more detail in subsection C. 



B. Off-Diagonal electron-phonon term: 
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FIG. 5. (a) MC snapshot of an 8 x 8 lattice at (n}=0.8 
A = and 7 = 0; (b) same as (a) but for A = 1 with mode 
Q^ 2 '; (c) same as (b) for A = 2; (d) Magnetic structure factor 
for the parameters in (b). 





FIG. 6. (a) MC snapshot of an 8 x 8 lattice at (n}=0.8, for 
mode Q , A = 1 and 7 = 0, after 2600 measuring sweeps; 
(b) same as (a) but after 3750 measuring sweeps; (c) same 
as (a) but after 4250 measuring sweeps; (d) same as (a) but 
after 5000 measuring sweeps. 





FIG. 7. Study of the effect of off-diagonal couplings, (a) 
MC snapshot of an 8 x 8 lattice at (n)=0.75, A = and 
7 = 0; (b) same as (a) but for 7 = 0.1 and mode Q' 2 ^; (c) 
same as (b) for 7 = 0.2; (d) same as (b) but for 7 = 0.6 . 
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FIG. 8. (a) Charge structure factor for different values of 
7 (strength of the off-diagonal EPI) and for A = 0ona8x8 
lattice at (n)=0.75, for mode (b) the magnetic structure 

factor for the same parameters as in (a) ; (c) the Drude weight 
for the same parameters as in (a); (d) the density of states 
for the same parameters as in (a). 

In this first exploratory study of the effects of the off- 
diagonal terms due to the EPI we will allow the parame- 
ter 7 in Eq.(8) and Eq.(10) to vary in the interval (0, 0.6), 
while gj (see Eq.(10)) will be kept equal to 1. In gen- 
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eral, we have observed that the effect of the off-diagonal 
term is to destabilize the stripes, since they become more 
dynamic. Examples of this effect can be seen in the snap- 
shots presented in Fig. 7 for (n) = 0.75 and 7 = 0, 0.1, 0.2 
and 0.6. The stripes become distorted as 7 increases and, 
eventually, AF domains separated by irregularly shaped 
hole-rich regions start to develop. 

In Fig. 8a, it is shown how the sharp maximum in the 
charge structure factor loses intensity as the stripes be- 
come more dynamic. However, notice that the magnetic 
structure factor (Fig. 8b) still shows incommensurability 
at q = (71", 37r/4) which is in agreement with neutron scat- 
tering data for dynamic stripes. For 7 > 0.4 a peak in 
5(q) develops at q = (0, 0). This seems to occur because 
the hole domains become ferromagnetic as it can be seen 
in Fig. 7c and d. 

The system also becomes more metallic since the 
Drude weight increases with 7 (Fig. 8c) and more spectral 
weight appears in the pseudogap in the DOS (Fig.8d). 

C. Diagonal and Off-Diagonal terms: 

When both diagonal and non-diagonal electron- 
phonon couplings are active simultaneously, we have ob- 
served that for values of A < 2, there is a competition 
between the localizing effect of the diagonal electron- 
phonon coupling A and the disordering tendency of the 
off-diagonal parameter 7. As a result, larger values of 7 
than in the case of A = are needed to destabilize the 
stripes. The stripe states are replaced by AF "clusters" 
separated by hole-rich regions as shown in Fig. 9. 




FIG. 9. (a) MC snapshot of an 8 x 8 lattice at (n) =0.875, 
for mode Q^ 2 \ A = 1 and 7 = 0; (b) same as (a) but for 
7 = 0.1; (c) same as (a) but for 7 = 0.2; (d) same as (a) but 
for 7 = 0.4. 



We also have observed that for A > 3 the modes 



and induce charge density wave (CDW) states. 
CDW domains are formed in order to accommodate the 
extra holes away from half-filling as shown in Fig. 10a 
for (n) — 0.875 and mode Q^ 2 \ This is the only case 
in which we have observed long range order developing 
in the lattice degrees of freedom. The dashed lines in 
the figure indicate different CDW domains. In Fig. 10b 
it can be seen how the off-diagonal coupling destabilizes 
the CDW state and a disordered state with bipolarons 
(two electrons trapped at the same site) is observed. 

As 7 increases the off-diagonal term opposes the trend 
towards bipolaron formation caused by the diagonal 
electron-phonon coupling and stripe-like structures reap- 
pear for some dopings. 




FIG. 10. (a) MC snapshot of an 8 x 8 lattice at (n) =0.875, 
for mode Q^ 2 \ A = 6 and 7 = 0. The dashed lines separate 
CDW domains; (b) same as (a) but for 7 = 0.2. 



The half-breathing modes, on the other hand, are not 
able to stabilize a CDW state even for large values of di- 
agonal couplings and they only produce disordered states 
with bipolarons when the diagonal electron-phonons cou- 
pling is large. Since CDW and superconducting states 
normally compete with each other, the half-breathing 
mode may enhance the development of superconductiv- 
ity. Unfortunately, off-diagonal long-range order which 
characterizes S or D wave superconductivity, cannot de- 
velop with adiabatic phonons. [34] Thus, this possibility 
should be explored for finite values of the frequency of 
the lattice vibrations in a full quantum calculation. 

D. Spectral Functions: 

In this subsection the properties of the spectral func- 
tions A(q, u>) will be discussed in detail. 

Numerical studies of the spin-fcrmion model without 
phonons in Ref. [21] have shown that the underdoped 
regime is characterized by a depletion of spectral weight 
along the diagonal direction in momentum space, where 
a very weak Fermi surface (FS) may exist. On the other 
hand, strong spectral weight, very flat bands, and a 
well defined Fermi surface (FS) are observed close to 
q = (7r,0) and (0, 7r). These results are in agreement 
with ARPES studies for LSCO, [35] material believed to 
have dynamic stripes, and also with numerical studies of 
models in which stripes have been built via a configura- 
tion dependent "stripe" potential in the t-J model [36] or 
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with stripe-like mean-field states. [37] ARPES results for 
LNSCO, where stripes have been observed, [27] indicate 
that the low-energy excitation near the expected d-wave 
node region is strongly suppressed. [38] 

In the absence of EPI the above mentioned char- 
acteristics are well reproduced by the SF model. In 
Fig. 11a the spectral function A(q, u>) along the path 
q = (0,0) - (tt,0)- ( 7r , 7r ) for ( n ) = 0.80 on a 12 x 12 
lattice in the absence of EPI is presented. A well de- 
fined quasi-particle peak is observed at (0, 0), a flat band 
appears close to (w, 0) and the peak crosses the Fermi 
energy at (iy,5tt/6) defining a hole-like FS. The spectral 
function along the diagonal of the Brillouin zone is pre- 
sented in Fig. lib. It can be seen how the spectral weight 
becomes very incoherent as the Fermi energy is reached. 
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FIG. 11. (a) Spectral functions along the path q = (0,0) — 
(7r,0) — (7r,7r) on a 12 x 12 lattice at (n}=0.8, A = and 
7 = 0; (b) same as (a) but along the diagonal direction in the 
Brillouin zone; (c) same as (a) but for A = 2 and mode Q^; 
(d) same as (b) but for A = 2. 
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FIG. 12. (a) Spectral functions along the path q = (0,0) — 
(7r,0) — (7r,7r) on a 12 x 12 lattice at (n}=0.8, A = 2 and 
7 = 0.4 for the half-breathing mode; (b) same as (a) but 
along the diagonal direction in the Brillouin zone. 

As a general rule we have observed that electron- 



phonon interactions increase the decoherence of the spec- 
tral functions, particularly close to the Fermi energy. Be- 
low, we will discuss the effects of the different phonon 
modes for the most interesting case of dynamic stripes 
(see, for example, Fig. 6). 

The diagonal electron-phonon coupling renders the 
system insulator, as discussed in subsection A. In Fig. 11c 
the spectral weight along the path q = (0,0) — (7T,0) — 
(it, 7t) is displayed for A = 2 and the half-breathing mode 
Q( 3 ) in a 12 x 12 lattice. It can be seen that the quasipar- 
ticles are less well defined and that a clear gap has opened 
at the Fermi energy. Along the diagonal of the Brillouin 
zone (Fig. lid) the spectral weight becomes more inco- 
herent and a gap also is observed. 

The off-diagonal EPI restores spectral weight in the 
gap but there are no quasiparticles close to the Fermi en- 
ergy as it can be seen in Fig. 12a and b where the spectral 
functions for A = 2, 7 = 0.4 and mode Q^ 3 ' are shown 
along q = (0, 0) — (n 7 0) — (it, n) and along the diago- 
nal of the Brillouin zone. For 7 finite and A = the 
only observed changes is that the peaks in Figs. 11a and 
b become broader. 

For A = 2 the breathing mode creates a larger insu- 
lating gap than the half-breathing mode and produces a 
larger decoherence of the quasiparticle peaks. The results 
for the shear mode are qualitatively very similar. 

The off-diagonal EPI for the shear and breathing 
modes tends to close the gap and increase the decoher- 
ence. 



IV. INFLUENCE OF PHONONS ON 
HOMOGENEOUS STATES AND GENERATION 
OF STRIPES 

Up to this point we have considered the effects of EPI 
on ground states that already presented charge inhomo- 
geneity. However, it is important to study whether the 
electron-phonon interactions proposed in this work can 
themselves generate charge inhomogeneity, in particular 
stripes, in a previously homogeneous ground state. 

In order to address this issue, we studied the S-F model 
with J = 1.5 instead of J — 2, value which was used in 
the previous sections (all the other parameters are kept 
the same). For (n) = 0.75 the ground state has an ho- 
mogeneous charge distribution as it can be observed in 
the MC snapshot presented in Fig. 13a. Note that despite 
the charge homogeneity this state presents magnetic in- 
commensurability due to a spiral spin arrangement in the 
vertical direction. 

One of the main results of this paper is our observation 
that a strong diagonal coupling with the shear mode gen- 
erates two horizontal or vertical stripes. An example can 
be seen in the snapshot presented in Fig. 13b for A = 2. 
In this case the holes act as boundaries between undoped 
antifcrromagnetic states and the magnetic incommensu- 
rability arises from the inhomogeneous charge distribu- 
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tion. A 7r-shift among the AF domains is observed as 
well. [39] 

The breathing mode also induces charge inhomogene- 
ity for a diagonal coupling A = 2. From the discussion 
in Section III-A diagonal stripes would be expected but 
we have observed two stripes with a zig-zag shape, i.e, 
the holes align diagonally at short distance scales but, 
on the whole, the stripe is still horizontal or vertical (see 
Fig. 13c). 




FIG. 13. (a) MC snapshot of an 8 x 8 lattice at (n)=0.75, 
J = 1.5, A = and 7 = 0; (b) same as (a) but for A = 2 with 
mode Q (2) ; (c) same as (b) with mode Q^ 1 ' . The dashed line 
indicates the stripes. 



This result indicates that diagonal EPI are able to in- 
duce stripe-like charge inhomogeneities in otherwise ho- 
mogeneous states. However, the S-F model also shows 
that although diagonal EPI stabilize stripes they are not 
necessary to induce them. Charge inhomogeneities can 
result even in the absence of EPI, just from magnetic 
interactions. 



V. CONCLUSIONS 

Summarizing, we have studied the effects of diagonal 
and off-diagonal electron-phonon interactions using non- 
biased numerical techniques on a model that has both 
striped and homogeneous ground states in the absence of 
EPI. We found that diagonal EPI tend to either gener- 
ate or further stabilize the existing stripes and turn the 
system into an insulator. Horizontal and vertical stripes 
are stabilized by half-breathing and shear modes, but the 
stripes generated by shear modes are more localized. The 
breathing mode stabilizes static diagonal stripes for in- 
termediate diagonal couplings. Large diagonal couplings 
stabilize CDW states for breathing and shear modes but 
this state is not observed with the half-breathing mode. 



On the other hand, off-diagonal electron-phonon cou- 
plings destabilize the stripes making the ground state 
more metallic, although non homogeneous. Instead of 
static stripes other kinds of inhomogeneous structures 
characterized by antiferromagnetic domains separated by 
barriers of holes are observed. 

The electron-phonon interactions undermine the quasi- 
particle peaks in the spectral functions close to the 
Fermi energy producing incoherent weight. Breathing 
and shear modes tend to open large gaps at the Fermi 
energy creating insulating behavior. The half-breathing 
mode, on the other hand, opens smaller gaps that are 
closed by relatively modest off-diagonal couplings, allow- 
ing for metallic behavior. 

Our results indicate that the half-breathing mode is 
most likely to play a role in non-insulating materials 
with vertical and/or horizontal "dynamic" stripes, such 
as LSCO, while the breathing mode should dominate on 
insulators with diagonal stripes like the nickelates. 

The most important result of this study is the fact that 
different phonon modes promote a diverse array of charge 
structures and that the relative strength of diagonal and 
off-diagonal couplings influences the transport proper- 
ties. Diagonal electron-phonon couplings promote insu- 
lating behavior, while off-diagonal interactions are crucial 
to achieve metalicity. It appears that the half-breathing 
mode off-diagonally coupled to the electrons is the most 
likely to produce non-insulating states with dynamical 
stripes as observed in the cuprates. This is in agreement 
with the experimental data which indicate the preva- 
lence of half-breathing modes in the high Tc cuprates. 
[6] The electron-phonon interaction introduces decoher- 
ence of the quasi-particle peak in the spectral function in 
agreement with ARPES measurements in LSCO. 

The crucial issue that remains to be explored is 
whether electron-phonon interactions are needed, in ad- 
dition to the magnetic exchange, in order to develop long 
range D-wave pairing correlations. Since superconductiv- 
ity arises from off-diagonal long range order, it cannot be 
generated with adiabatic phonons. The next step will be 
to study off-diagonally coupled half-breathing modes at 
finite frequencies. 
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